Chiral Symmetry Breaking in Bent- Core Liquid Crystals 
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Abstract 

By molecular modeling we demonstrate that the nematic long-range order discovered in bent- 
core liquid crystal systems should reveal further spatially homogeneous phases. Two of them are 
identified as a tetrahedratic nematic (Nt) phase with symmetry and a chiral tetrahedratic 
nematic (N T ) phase with D2 symmetry. These new phases were found for a lattice model with 
quadrupolar and octupolar anisotropic interactions using Mean Field theory and Monte Carlo 
simulations. The phase diagrams exhibit tetrahedratic (T), Nt and N T phases, in addition to 
ordinary isotropic (I), uniaxial nematic (Njj) and biaxial nematic (Nr) phases. To our knowledge, 
this is the first molecular model with spontaneous chiral symmetry breaking in non-layered systems. 
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The chirality of liquid crystals is commonly attributed to the presence of optically active 
(chiral) molecules Such intrinsic molecular chirality distinguishes right-handed and left- 
handed molecules leading to a rich variety of chiral phases, where cholesteric phase [l] and 
blue phases j^j serve as typical examples. 



Interestingly, this classical liquid crystalline chirality s 



core (banana-shaped) liquid crystalline molecules |3|, 4], |5|, |6|, |7|, |8|, |9_|, [10|, lU| • These fascinating 
compounds have substantially different physical properties than the molecules of the tradi- 
tional calamitic materials, and have recently been a source of discoveries that are changing 
our view on 

Several observations indicate that achiral bent-core molecules can acquire conformational 
chirality in layered (smectic) mesophases 0, [k]] . They can form ferroelectric and homoge- 



rould be revised for achiral bent- 
9. 



neously chiral states (4], or coexisting left and right-handed chiral domains in a nematic phase 
6|, |8j. All basic theoretical interpretations of these observations introduce at a macroscopic 
scale a tetrahedratic order parameter field. The interplay between this order parameter and 
the (nematic) quadrupolar tensor order parameter 12|, ll3[ appears to be essential for the 
emergence of the chirality of the achiral bent-core molecules |7|, |9|. In the bent-core molecu- 
lar systems this coupling generates not only the most interesting chiral symmetry breaking. 
The new tetrahedratic (T) and tetrahedratic nematic (Nt) phases should also be stabilized 
Recently, the experimental indication of stable T and Nr_ phases was indeed reported 



14j , along with earlier discovery of the biaxial nematic phase |15l . 

This letter analyzes some consequences of molecular interactions that have tetrahedratic, 
uniaxial and biaxial components. For a complete understanding of the macroscopic proper- 
ties of mesophases in the bent-core systems, a first-rank tensor generating structures with 
polar order, is also required. To simplify the analysis, the microscopic model being discussed 
is limited to quadrupolar and octupolar interactions. In addition to isotropic and nematic 
phases, the model is shown to stabilize T, Nt and chiral tetrahedratic nematic (N T ) phases. 
These results are what we believe a first microscopic demonstration of spontaneous chiral 
symmetry breaking in non-layered bent-core liquid crystals. 

We assume the molecules to occupy the sites of a three dimensional simple cubic lattice 
with interactions limited to nearest neighbors (coordination number: d — 6). The total 
interaction potential is then given by H = | X^<t • > V(fij, fij), where V(Qi,Qj) is the 
orientational interaction between a pair of molecules (or molecular complexes in general) 
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and < i,j > denotes the nearest-neighbor molecules i and j. Our calculations are based 
on the simplest attractive, (9(3)-invariant interaction V(fij, ttj) that involves second-rank 



quadrupolar tensor with uniaxial (T^) and biaxial (T^) components, and a third-rank 
octupolar tensor (T2 ). The tensors are built out of the orthonormal tripod Q k of vectors 
{a. k ,b k ,c k } defining the orientational degrees of freedom of the k-th molecule. Under these 
assumptions V(Qi, Slj) takes the form 
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and where T $ ■ T { 3 



17|| . The summation in Tj, 3 '* runs over all permutations of 



(a&, bjfc, Cfc). The symbol ' • ' denotes the scalar product formed by a full contraction of 
the Cartesian indices, and T^ffij 
rotation matrices 



^mm' are linear combinations of Wigner 



The interaction (1) can be interpreted within point dispersion forces approximation. The 
quadrupolar tensor + v^AT^ is then proportional to the anisotropic part of the di- 
electric polarizability tensor of a molecule while T/j, is the T^-symmetric component of the 
third-order polarizability tensor. The model parameters {e, A, r} are related to the molecular 
absorption frequencies. 

Special cases of the model have already been studied. For r = A = the potential in 
Eq. (1) reduces to the well-known Maier-Saupe or Lebwohl-Lasher [jjj] potential, which 
describes a phase diagram with isotropic and uniaxial nematic phases connected by a first- 
order phase transition. Second is the one where A ?z and r = 0. Potential (1) then 



reduces to the model proposed by Luckhurst et al. (20|, which was extensively studied by 
Biscarini et. al. 21[. The model describes a phase diagram with uniaxial nematic and 
biaxial nematic phases connected by the second-order phase transition. The phases include: 
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a prolate uniaxial (Nu + ) phase, an oblate uniaxial {Nu_) phase, a biaxial (Nb) phase, and 



an isotropic (I) phase. A self-dual point for which A = 1/ 18|, |2l( separates a phase in 
which the molecules are of distorted prolate form (A < l/v6) from a phase in which the 
molecules are of distorted oblate form (A > l/y/E). Bates and Luckhurst have proposed 
a simple relation between e, A and the opening angle for bent-core molecules based on 
segmental second rank interactions 2^ ]. 

When only the term proportional to er is retained in Eq. (1), the resulting model cor- 
responds to a purely tetrahedratic coupling. It was introduced by Fe l |23l| and studied via 
Mean Field (MF) and Monte Carlo (MC) simulations by Romano [24]]. For this model, 
MF theory predicts a second order phase transition from I to T phase and MC simulations 
indicate a weak first-order transition from I to T. 

Combining quadrupolar and octupolar interactions in Eq. (1), yields new possibilities 
of which the most notable one is the spontaneous breaking of chiral symmetry. Such a 
symmetry breaking is already manifested in the ground state properties of the interaction 
of Eq.(l). Indeed, to be consistent with Eq.(l), the average molecular configurations of two 
bent-core molecules in a chiral phase must belong to one of two configurations of opposite 
chirality (0 < 5 < vr/2), as shown in Fig. 1. Due to the global 0(3)-invariance of H, the 
two configurations of different chirality are the source of two homo-chiral domains, which 
are present with equal probability in the system's configuration space. When r = (8 — 0) 
or A = (8 = 7r/2) the molecular configurations are nonchiral and can be brought into 
coincidence by a rotation. In particular, these configurations produce Njj, Nb, T and Nt 
phases. In order for both chiral configurations to be equivalent, the degrees of freedom of 
the k-th molecule should involve a rotation and parity pk = B-k ■ ^b^. x c^j = ±1, thereby 
reflecting an 0(3) symmetry. Consequently, the free energy for the system composed of N 
such molecules is given by (3F = — InZ, where Z = Ylk=i (§Sp fc =±i / d^/U ex P [ — flH], 
and (3 = l/iksT) is the inverse temperature. 

We will apply two methods to determine phase diagrams from F. The methods are 
the MF approximation and the Metropolis MC simulations with a crucial step being the 
identification of the order parameters. The identification is achieved by expanding the one- 
particle distribution function P(p, Q) in a series of symmetry adapted, real A-functions: 
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P (& n )= E ^Ps^n'Ps^UV), (3) 

L,m t m' ,s 



with x = \ j: p=±1 j x(p, n)P(p, n)dn, Pa = (i, P ), and 

5 E / dO^A^(^) Ps>A%? n ,(n) = ^^5 LL ,5 mml 5 nnl 5 ss , (4) 

For the interaction of Eq. (1), the allowed symmetry reduction is shown in Fig. 2 along 
with the primary order parameters that acquire non-zero averages when crossing various 
phase transitions. In order to be consistent with L = 3 primary tetrahedratic order param- 
eter, we list all remaining (secondary) molecular order parameters of rank L < 3. These 

(2) (2) (3) (2) (2) 

are: (a) A 02 for Ny; (b) A 20 and pA 22 for Nb', (c) pA 2( j an d pA 22 > a l° n g with Njj secondary 
order parameters, for Nt] and finally (d) for N T all aforementioned order parameters are 
nonzero, in addition to p, pA^ , and pA^ 2 . Here 

/ i \2+5om+5 0m i 

where X> are the Wigner rotation matrices. 

We are now in position to evaluate the MF approximation to F and obtain the equilibrium 
properties of our system. We have identified six stable phases of the model (1) that cover 
the whole symmetry reduction flowchart of Fig. 2. The phases are: (a) the isotropic phase; 
(b) the uniaxial prolate or oblate nematic phase; (c) the biaxial phase; (d) the tetrahedratic 
phase; (e) the prolate or oblate tetrahedratic nematic phase and (f ) the chiral tetrahedratic 
nematic phase. Exemplary phase diagrams for r = 1 and and for r = 9 are shown in Figs. 3 
and 4, respectively. Fig. 3 shows a phase diagram in which the high-temperature region is 
dominated by the nematic phases. Fig. 4 shows a phase diagram in which occurs a direct 
I — T phase transition. Phase diagrams for intermediate values of r can partly be deduced 
by extrapolation. One of them is of particular interest. Namely, for r = || six phases: /, T, 
Nt+, Nt*, Nt- and Njj_ meet at a single multicritical Landau point. All phase transitions 
found are of second order except when more than one order parameter acquires a nonzero 
value at a bifurcation. Such first order phase transitions occur at / — Nt, T — Nt or I — Njj. 
A fuller account of the properties of this model is deferred to our future publication. 

To check the validity of the MF predictions the phase diagrams were also determined 
from MC simulations. The simulations were performed on a 16x16x16 lattice with periodic 
boundary conditions. Each MC move included a rotation of a molecule's orientation and a 
parity inversion. The size of MC rotational moves was selected to produce an acceptance 
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ratio between 30% and 40% in the ordered phases. Typically, 50 000 to 200 000 lattice sweeps 
were used to thermalize the system and 60 000 to 200 000 sweeps for measurements. In the 
MC simulations, phase transitions were detected by observing the temperature dependence 
of the order parameters. In turn, these quantities were determined from the asymptotic 
behaviors of the correlation functions G^ m ,(\i— j\) = (f2$) • T^/(flj) and G pp (\i—j\) = 
pip] for large \i — j\. Here the overline indicates an ensemble average. In the simulations, 
the temperature resolution 5t satisfied 5t = 0.01. At the above resolution (and system's 
size) our simulations were unable to distinguish between weakly first-order and second-order 
phase transitions. The MC transition curves lie below those obtained from MF. The MF 
results worsen from a discrepancy of 8% to about 20% as A increases from A = (/ — Njj) 
to A = 0.3 (Nt — N T ). The I — Nu part of our MC diagram agrees with MC results of 



Biscarini et al. 



2l|. 



The new liquid crystal phases appear due to the combination of octupolar interactions 
with quadrupolar contributions. In addition to a tetrahedratic nematic phase, existing in 
prolate (Nt+) and oblate (Nt-) versions, the N T phase appears stable. It does not appear for 



A = 0orA = yf, which correspond to a uniaxial limit. From the leading MF contribution 
to the parity order parameter, p, an induced homogeneous chirality is seen to emerge when 
both biaxial and tetrahedratic orders condense. The estimation gives 

V2 



V 



2 10t 4 



r\ (-3 + 2A 2 ) Ag (Ag + v^AAf) , (6) 



which vanishes when biaxial or tetrahedratic order parameters are zero. A maximal transi- 
tion temperature to a chiral phase is achieved for the self-dual point of A = 1/ \/6- 

To conclude, bent-core liquid crystals can stabilize the elusive thermotropic biaxial ne- 



matic phase 
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and reveal a path to a series of new spatially homogeneous but 
anisotropic liquids |7|]. Of these liquids the most interesting one is the N T phase, which 
should emerge from a nonchiral liquid as a result of spontaneous chiral symmetry breaking. 
The simplest 'spin-spin'-type of molecular model of Eq. (1), with quadrupolar and octupo- 
lar interactions supports this scenario, but, unlike chiral phases of ordinary chiral materials, 
the orientational order in N T appears spatially uniform. It is worth noting, however, that 
inclusion of higher-order cross-coupling terms between uniaxial, biaxial and tetrahedratic 
interactions can superimpose a spatial modulation to Nt and N T [17J. Thus the struc- 
tures stabilized in our model can serve as the long- wavelength limits to a family of spatially 
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modulated chiral structures that can possibly condense in the presence of higher-order in- 
teractions. Owing to general form of the interaction, Eq. (1), these conclusions should apply 
to any system where tetrahedratic and quadrupolar order may simultaneously coexist. 
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FIG. 1: (Color online) Ground state configurations of opposite chirality (0 < 5 < vr/2) for a pair 
of bent-core molecules in Nj, phase; the configurations for < 5 < tt/2 cannot be brought into 
coincidence by a proper rotation. (a m ,b m ,c m ) are mirror images of (a, b,c). 
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Isotropic: I { 0(3) } 



Uniaxial: 



Biaxial: 

N B lD 2h ) 




Tetrahedratic: 



oo 

Nematic 
Tetrahedratic: 

N T { D 2d } 



Chiral: N T * { D 2 } 



FIG. 2: (Color online) A flowchart of phase transitions between liquid-crystal phases for molecules 
interacting through pair potential, Eq. (1). Primary order parameters, which become nonzero at 
the transitions, their symmetry groups and abbreviated notation for the structures are indicated. 
In addition, sketched are one-particle spherical distribution functions P ~ const + al20 + ^(^22 + 
Y2-2) + ci(Y^2 — Y3-2) illustrating symmetries considered. 
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FIG. 4: (Color online) MF phase diagram for r = 9 (see captions to Fig. 3). Simulations here 
are difficult since parity degrees of freedom condense to a glassy state. The Nj, phase cannot be 
reached by standard Metropolis simulations. 
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